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ABSTRACT 


Hydrodynami cs of an axi -“syimnetri c combination blowi 
steel making process has been analysed mathematically. Tbis 
reci rculatory flow field within the reactor vessel has beer 
represented interms of the turbulent Navi er -Stokes equations ir 
conjunction with the k-^ turbulence model and the resultant set of 
partial differential equations together with the appropriate 
boundary conditions solved via a control volume based finite 
difference procedure embodied in the popular TEACH-T computer 
code. Prior to any evaluation of the model predictions » sens! tivity 
of computations to the choice of residual value » noda 
configurations and differencing schemes C e. g. > upwi nd , hybrid, 
exponent! al , power law etc. have been assessed and results 
independent of these established. 

The bottom/submerged gas injection phenomena in a 
combined blowing reactor has been represented via the single phase 
variable density formulation technique, wherein the rising two 
phase gas-liquid mixture has been treated like a homogeneous fluid 
of reduced density. The gas volume fraction and hence the density 
of the two phase mixture have been estimated from a knowledge of 
the dimensions of the two phase region and the latter average I 
rise velocity. To this end, several calculation procedure^ 
Cviz. , formulations based on slip vs no-slip, constant gas volume! 
fraction vs variable gas volume fraction in the gas -liquid plume 
etc. 3 reported in the literature have been critically examined. 
Through extensive comparisons between theory and experiments it is | 
demonstrated that the bottom injection phenomena can be ! 
mathematically adequately represented by assuming gas/bubble 
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slippage and considering a constant rise velocity of the two phase 
mi xtur e. 

In representing the impinging gas jet Ce. g. ,the top gas 
injection system!) the shearing force imparted by the outwardly 
flowing gas across the cavity surface Ce. g. ,the gas -J i quid 
interface!) has so far been assumed to be the only mechanism 

producing the bulk liquid recirculation. However » the present 
investigation appears to suggest that flows generated through such 
considerations alone are not sufficient enough to drive the flow 
recirculation observed experimentally. Consequent! y* as an 
additional mechanism to the shearing force* the force of 

impingement estimated from the Increased potential energy of the 
melt Cviz. * deduced from the corresponding increase in bath 
height !> has been considered and the top blowing system 
represented accordingly via an approach based on these two 

mechanisms. Flow fields thus predicted have been evaluated 

against equivalent experimental measurements on water models and 
excellent agreement achieved. Finally* hydrodynamic models of top 
and bottom injection systems have been coupled and a steady state 
fluid flow model for combination blown steel making system 
proposed. Numerical predictions from the model have also been 
evaluated against experimental measurements and reasonably good 
agreement between theory and experiments 


demonstrated. 



Nomenclature 


iii 



Coefficients of discretization equation 

representing the effect of convection and 
diffusion 

Summation of all neighbour point 

coef f i ci ents 


Constant defined in CI . 43 : has a 

value of 7 

Dissipation rate constant : has a value 
of O. 09 

Constant defined in CI.53: has a value 
of O. la 

Constant defined in CI . 43 : has a value 
of 0.107 

Constant defined in CI.63: has a value 
of 9. 1 

Constant defined in C I 1.103: has a value 
of 1 . 43 

Constant defined in C I 1.103: has a value 
of 1 . 93 

Diameter of a lance nozzle » m 
Diameter of a gas bubble* m 

~S 

Acceleration due to gravity* m s 
Height of the liquid bath* m 
Distance between analytic origin and the 
base of the vessel , m 



Constant defined in Cl . 4Z> : has a value 
of 1. O 

Turbulence kinetic energy, m s 

Distance between orifice of lance nozzle 

and free surface, m 

Pressure within the liquid, Pa 

3 

Gas flow rate, m /s 

the radial coordinate, m 

Radius of the bath, m 

Radius of the two phase plume, xn 

Sum of residual source 

Constant plume velocity, m/s 

Bubble slip velocity, m/s 

Mixture rise velocity in the two phase 
r egi on , m/s 

Axial velocity component, mys 
Radial velocity component, m/s 

Volume fraction of the bulk phase in a 
c on t r ol vol ume 

Volume fraction of the two phase region in 
c on t r ol vol ume 
the axial coordinate, m 
Gas void fraction 

Gas void fraction at the axis of symmetry 
Shear stress acting across the cavity 
surface, kg m ^ 



V 



a 

e 


0 

X 


Density of gas » kg m 

~3 

Density of liquid, kg m 

-1 -1 

Viscosity of the gas, kg m s 

“1 -1 

Effective viscosity, kg m s 

Turbulence viscosity, kg m s 

~1 -1 

Laminar viscosity, kg m s 

Turbulence kinetic energy dissipation rate 
per uni t mass , m s 

Constant in the ks model : has a value of 

1. O 

Constant in the k-^: model : has a value of 

1 . 3 

Any general variable 

Half jet cone angle, radians 



List Or Figures 


Vi 


Figures 
1 . 1 

II . 1 


II . B 


II . 3 


II . 4 


II . 5 


II . 6 


III . 1 


The combination blown steel making system. 

Schematic diagram of the idealised 
combination blown s>'^tem used for the present 
mathemati cal r epr esentaii on , 

The continuum approach used to define the 
density distribution in the calculation domain. 
Segmentised cavity contour and the 
procedure adopted to evaluate components of 
shear force/stress at any point on the 
cavity surface. 

Representation of the various boundary 
conditions applied to the numerical 
solution of the governing partial differential 
equations. 

Various t>'pes of control volumes and the 
staggered grid system applied to the present 
numerical solution scheme. 

Flow sheet of the TEACH-T computer 
program 

Sensitivity of computations to the choice 
of residual CNARSD values illustrating 
the radial variation of the vertical velocity 
component of flow for three typical NARS 

f=> o 

values CQ=1. 67x10 m /s, R=0. 15m, L=0. 21m 15x10 


grid s>'stenO. 




vii 

III. 2 

Sensitivity of computations to the choice 

of nodal configurations illustrating the radial 

variation of the vertical component of 

flow for three different grid systems 

C Q=1 . 67x1 , R=0. 1 5m , L=0. 21 m , NARS=5xl 0~^:» 

III . 3 

Sensitivity of numerically predicted 

vertical velocity component in a water model of 

27 

a gas stirred ladle illustrating the 

influence of differencing scheme employed on 

predi cted resul ts. 

III . 4 

Theoretically estimated centre line gas 

voi dages i n a water model 1 adl e and thei r 

comparison with equivalent experimental 

measurements . 

III . 5 

Theor etl cal 1 y esti mated cr oss-secti onal 

average gas voi dages in a water model and their 

comparison wl th equi val ent exper i mental 

measur ement s . 

111.6 

Theoretically predicted vertical velocity 

component within the rising plume and their 

comparison with experimental measurements 

of Lehner et al . ^ ^ . 

III. 7 

The shearing forces acting at a location 

on the cavity surface and the configuration 

of the shear layer tick ness with 

reference to the finite difference grid system. 


v/ll 


III . 8 


III . 9 


III . lO 


III . 11 


1 JI . IB 


III . 13 


Numerically predicted vertical velocity 

component at different depths and their 

comparison with experimental measurements of 

31 

Jagannaihan et al . 

Sensitivity of numerical predictions to 
the values of 1^ illustrating the 
variation of vertical velocity component 
at two different depths in the model reactor. 
Sensitivity of numerical predictions to 
the choice of shear layer thickness, illustrating 
the variation of vertical velocity 
component at two different depths in the 
model reactor. 

Numerically predicted \'ertical velocity 

component at different depths and their 

comparison with experimental measurements of 

31 

Jagannathan et al . 

21 

Experi mental 1 y measur ed vel oci ty f i el d 

in a water model of a combination blown 
steel making system. 

Numerically predicted flow field in a 

water model of a combination blown 

reactor C conditions corresponding to 

21 

those of Zhang and coworkers 3. 


IX 


III. 14 The configuration of the cavity surface 

with reference to the finite difference grid 
illustrating the principles involved in 
the cell porosity method applied to 
exclude the cavity volume from the domain 
of calculation. 

III. 15 Theoretically predicted vertical velocity 

component at two different depths in the 
B1 

model \'essel illustrating the effect of 

consi deri ng/excl udi ng the cavity volume 

on the computed results. 

Predicted turbulence kinetic energy 
2 ~2 

Cm s D distribution in a water model of a 
combination blown reactor C corresponding 
to the experimental conditions of Zhang and 
coworkers 3. 


Ill . 16 


X 


List Of Tables 


Tabl e 1 1 . 1 
Tabl e I I I . 1 

TabJe III, 2 


Table III. 3 


Tabl e I I I . 4 


Values of constants in the turbulence model 

Experimental conditions of different 
investigators used in the computation of 
results presented in Figs. 1 1 1. 1 through III. 5. 
Theoretically predicted mean speed of 
liquid recircul ati on via various 

differencing schemes and their comparison 
wi th expert mental measur ements of 

Mazumdar et al . . 

Experimentally measured mean speed of 

C273 

liquid recirculation and their 

comparison with those estimated 

theoretically to illustrate the sensitivity 
of prediction to the choice of gas volume 
fraction cal ul ati on procedures. 

Experimental conditions of different 
investigators used in the computation of 
results presented in Figs. III. 7 through 


III. 16. 


Chapter I 


INTRODUCTION 

1.1 Introduction To The Thesis 

The injection of reactive and non~reaciive gases into 
melts contained in metallurgical furnaces, ladles etc. is a 
frequent phenomena encountered in numerous metal processing 
operations. In either case, the injected gas promotes fluid motion 
within the vessel and thereby exerts considerable influences on 
phenomena, such as, slag-metal react! ons , thermal and particulate 
homogenisation, inclusion float out or re-entrainment etc.. 
Consequently, gas injection induced stirring in metallurgical 
reactors has since long been the subject of considerable interest. 
Moreover , thi s has been recognized to be of vital importance in 
determining the efficiency of numerous associated processes 
Cviz. ,melting/dissolution, dispersi on/mlxi ng etc.!> carried out in 
gas stirred reactors. 

Typically, two different modes of gas injection have 
been applied to produce the required stirring and these include ; 
Ci3 gas Jet impinging on the surface of metallic melts 
and 

CiiDsubsurf ace injection of gas through a lance, tuyere or 

porous -pi ug. 

As an example to the first category, mention can be made of the 
top bl own oxygen steel making systemC vi z . , the LD pr ocessD whi 1 e 
for the latter, depending on the intensity of the blow a wide 
range of metal processing operations Cviz. , from bottom blown 
steelmaking processes to ladle processing of steelD may be 
cited. It has been realised for some time that the intensity of 
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bath stirring In top blown reactors Is relatively sluggish, as the 
transfer of momentum from the impinging gas stream to the bulk 
liquid phase is rather Inefficient. On the contrary, bottom 
injection provides vigorous bath agitation and thus exacerbates 
the rates of various mass transfer controlled phenomena, 
characteristic of high temperature operations. Indeed, in 

combination blown steel making processes desired metallurgical 
results have been successfully obtained by judiciously combining 
the top C this supplies the necessary oxygen for metalloid 
oxidation react! onsD and the bottom/submerged gas injection 
si'stems Csee Fig. I.l>. Since the primary objective of the bottom 
gas injection s>'stem in a combined blown reactor is to provide 
necessary stirring rather than supply oxygen C the latter is known 
to be detrimental, since the bottom/tuyere life decreases Ly 

injecting oxygen through the bottom, affecting considerably the 
turnover rates per converter^ , consequently injection of inert 
gases like argon etc. have been typically practised. Furthermore, 
the "gas-liquid coupling" in submerged gas injection being more 
efficient Cthan an impinging gas jet systenD , requires a small 
volume of gas to be injected in relation to the top blown gas 
volume C typically about 5?0 to produce the required stirring 
ef f ect . 

Since the flow fields induced in combined blown reactors 
Csay, the play paramount roles in controlling the 

associated processes, such as slag-metal reaction rates, 
melting/dissolution of scrap, refractory lining wear etc. 
consequently, adequate knowledge of bath hydrodynamics is critical 
to the overall process anal>'sis. However, since high temperature, 
visual opacity and relatively large sizes of metallurgical 
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Refractory lined 
vessel wall 



Fig. 1.1. The combination blowing steel making 
system . 





react-ors pose serious limit, at. Ions to direct experimental 
observations, consequently, an alternative approach such as cold 
model stud! es/mathemati cal modelling C or both In conjunction D 
has to be adopted to develop understanding of fluid flow and 
other associated phenomena. 

In the present work, hydrod>Tiamics of a combination 
blown steel making process has been modelled mathematically and the 
model predictions assessed against experimental work reported in 
the literature on equivalent room temperature physical models 
Ce. g. , water modelsD. Towards these, at the onset, a brief review 
of relevant literature has been presented and the scope of the 
present work discussed. In subsequent chapters, the development of 
the mathematical model and the numerical modelling techniques have 
been outlined in detail and finally, several computed results have 
been compared against experimental observations. Summary and 
conclusions as well as recommendations for further work appear 
respectively in the concluding chapters. Computer programs etc. 
have also been assembled in an Appendix. 
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I#E Review Of Relevant Literature 

EXiring the last two decades or so numerous 

investigations have been carried out on the combination blown 
steel making process and phenomena such as » mixing, slag-metal 
reactions, mass transfer rates etc. have been extensively studied. 
However, fluid flow aspects of the combination blown process has 
remained relatively less well investigated. Furthermore, little or 
no effort has so far been made to address fluid d>'namics of such 
systems from mathematical points of view. Consequently, barring 
only one instance, no studies on mstthematical modelling of 
combination blown process have been reported in the literature. 

Similarly, mathematical model studies of the top gas injection 
system Cviz. , the LD configuration!) has also been rare and only 
one such study so far have been noted. In contrast, however, a 
large number of investigations on the mathematical modelling of 
submerged gas injection s>^tems appear in the literature. Since 
combination blown process embodies both top and bottom injection 
systems, consequently, fluid model studies on these systems have 
also been included in the present review in addition to those on 
combination blown system. 

The bottom/submerged gas injection system in the 

combined blowing reactor invol\^s the Injection of gas into the 

melt through one or a number of tuyeres located at the base of the 
vessel. Since the volume of gas injected through the bottom is 
typically small, the injection process consequently, leads to 
phenomena si ml la to those observed in ladle processing of steel. 
The Injected gas in the immediate vicinity of the tuyeres rapidly 



devolves into plumes of spherical cap bubbles* rise due io 
buoyancy* imparting energy to the bulk liquid and leaves the 
vessel through the free surface. Therefore, it is more appropriate 
to consider a buoyant plume rather than a Jet model to describe 
the bottom injection phenomena. Thie interactions between the gas 
and the surrounding liquid induces flow recirculation in the 
vessel as has been mentioned already and this facilitates many 
metallurgical processes, such as, slag-metal reactions, 
dissolution, mixing etc. . 

There are two different computational procedures 
typically been adopted to mathematically describe the 
hydrodynamics of such two phase buoyancy driven systems and these 
i ncl ude: 

CaD models derived on basis of Euleri an frame of 
reference* which are further classified into two groups viz. , 

CiD the two phase models 

CiiDthe quasi -single phase models 

and 

CbD models derived on the basi^ of Lagrangian frame of 

ref er ence. 

In the two phase modelling, two sets of continuity and 

r 1 

momentum equations are solved, one for each phase .This approach 
of modelling being relatively more complex from the computational 
view points, has therefore, received less attention in the 
literature. Momentum equations Tor individual phases are 
interlinked via the interface friction f actorCe. g. » slipD and in 
addition, no assumption is made about the gas volume 

fraction C vi z . , voi dageZ) distribution in the plume. The latter is 
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rAih©'r malted -through the solution of the gas phase continuity 

and momentum equations. 

The quasi -single phase modelling invol\^s solution of 

only one set of continuity and momentum equations. Such an 

approach* however* requires a prior specification Ceither 

estimated or empi ri cal Z)of the plume geometry Cviz. , the gas- liquid 

zoneD as well as the gas voidage. The buoyancy force per unit 

volume Cp^gcO estimated on the basis of the assigned gas voidage 

Cviz. *cO is added to the vertical direction momentum equation as 

the required driving force. Eulerian single phase modelling has 

received considerable attention in the literature because of its 

relative simplicity and ease of computational efforts. This 

t echni que , f or exampl e * has beer used extensi vel y by E>ebr oy et 

al . , Szekely and coworkers * Sahai and Guthrie and 

C 7 8D 

Mazumdar and Guthrie * to describe turbulent reci r cul ator y flow 
fields in bubble .stirred ladles. 

In the latter approach * C e. g. *CbDD the second 

phased , e. * the gasD is treated essentially by a Lagrangian 
technique* in which a bubble trajectory equation is used in 
addition to the liquid phase continuity and momentum equati onsC the 
latter are same as those in the quasi- single phase modellings. 
The bubble trajectory equation together with the momentum 

equations provide estimates of the two phase region, gas voidage 

^ C9S 
etc. 


Of the three approaches of modelling a submerged gas 
injection system mentioned above* the quasi -single phase modelling 
technique has been more popular. Consequently, the present 


1 i ter atur e r evi ew 


the mathematical modelling of the 


on 
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bolicjm -^submerged gas Injection system has been restricted to the 

studies based on such an approach only. 

CS3 

Etebroy, Majumdar and SpaJ di ng were the first to 

introduce the concept of a buoyancy driven plume in the 

metallurgical ladles. They considered the density deficit between 

the gas and the liquid phases to be the driving force Cviz. , the 

buoyancy^ producing the flow recirculation and modelled this 

effect through adding the buoyancy force term to the axial 

C 2D 

direction momentum equation. The authors adopted the Eulerian 

single phase modelling approach Ci.e. .one set of continuity and 

momentum equationsD together with a bulk effective viscosity 

model for describing the turbulent flow in the s>'stem. 

Symmetry conditions were applied to all the variables C si nee the 

configuration considered was axi -symmetricD . In order to estimate 

the density deficit Cor the void fraction^ in the two phase 

region, zero slip condition Ci.e., the gas and the liquid move 

together^ has been applied. Numerical predictions of flow 

variables were assessed against experimental measurements of 

CllD 

Szekely et al . and qualitative agreement between theory and 

experiments have been demonstrated. 

C3D 

El-kaddah and Szekely proposed a fluid flow based 

mathematical model for estimating sulphur partitioning between 

slag and the metal phases in a gas stirred ladle. The 

desulphurisation kinetics was described through prediction of 

turbulent flow fields via the solution of Navier-Stokes equation 

C 1 2D 

in conjunction with the k-s model. The calculation procedure 


in addition involved prior specification of gas voidageCcO through 
the solution of an ordinary differential equation, viz. , 



. . . Cl . i:) 


dz 


IE CtanX:) u z 3 
z 


and Incorporating in an expression, viz. , 


aCzZ> = 


tQ -nr aCl-cO U 3 
g c s 


En f r U dr 

o-' 


... Cl . ED 


in which is the slip velocity. Fluid flow based predictions 

were compared with the observed rate of desul phuri sati on deduced 

from high temperature melts and reasonably good agreement between 

theoretical and experimental results has been reported. 

Grevet * Szekel y and El -Kaddah subsequently reported an 

extensive fluid flow study in a water model of an argon stirred 

ladle. In this investigation, Eulerian single phase modelling 

together with an assumption of bubble slippage have been 

incorporated. The formulation of governing equations etc. were 

essentially identical to those reported earlier by El -Kaddah and 
C 33 

Szekely . Excellent agreement between predicted and experimental 

velocity fields C deduced through Laser Doppler VelocimeterD was 

achieved. However, predicted and experimental turbulence kinetic 

C43 

energy were only in moderate agreement. The authors 

acknowledged that such discrepancies were perhaps associated with 

the incorrect selection of constants or with the Inherent 

shortcomings of the k-£: turbulence model for representing such 

C 43 

reel rculatory flow systems. Furthermore, the authors suggested 
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that the principal mechanism of momentum transfer In bubble driven 
systems Is due to fluid convection^ C rather than turbulent 
diffusionD and thus predictions regarding the velocity fields were 
not expected to be very sensitive to the choice of a turbulence 
model . 

C 5 

Sahal and Guthrie ' through their extensive 

experimental and theoretical studies, postulated that bubbles 
forming at nozzles or orifices rapidly devolve into a plume of 
spherical cap bubbles and that this is true for many gas injection 

processes carried out in holding or transfer vessels. The rising 

bubbles Impart energy to the surrounding liquid and produce flow 
recirculation within the vessel. Equating the total energy 

supplied by the bubbles to the turbulence energy dissipation 
losses in the bath, a simple macroscopic expression for estimating 
the average rise velocity of the bubble plume has been proposed 
according to ; 


q1/3 ^1/4 


U = 4. 4 C D ... Cl . 3D 

P ^1X3 

In subsequent numerical prediction of flow etc. , they adopted the 

above relationship to estimate the gas voidage Cthey considered no 

slip conditionD in the two phase region. Furthermore, fixed 

velocity at the axis of the vessel has been applied as the 

required boundary condition instead of the typical zero gradient 

conditions. Predicted values were compared with experimental 

measurements of their own as well as with those of Lehner et 
Cl 3D 

al . and reasonable agreement between theory and experiments 


has been demonstrated. 
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C V 0Z) 

Mazumdar and Guthrie * proposed hydr od>Tiami c models 

for various axi symmetri c gas injection configurations found in the 

industrial ladle refining operations. The authors estimated gas 

void fraction via equation Cl . 3D . The hydrodynamic model in 

essence was based on the solution of turbulent Navi er -Stokes 

equation together with the k-s turbulence model. EXie to some 

discrepancy between theory and predictions for the CAS 

C4D 

configuration the empirical constants in the ks model were 

changed by some adhoc values. Numerically predicted values Cflow 
and turbulence parametersD thus obtained were shown to be in 

excellent agreement with experimental results. 

Cl 4D 

Salcudean and coworkers proposed a three dimensional 

calculation procedure to describe non-axis>'mmetric gas injection 
configurations through the solution of turbulent Navi er -Stokes 
equation, incorporating the k-s turbulence model. There were, 
however, large discrepancies between theoretical estimates and 
experimental measurements particularly in the immediate vicinity 
of the nozzle. Even though good agreement between predicted and 
experimental values were shown for elsewhere, the authors 
suggested that the discrepancy could be from numerical errors due 
to false diffusion phenomena Ce. g. , related to the use of the 
differencing scheme etc.D as well as due to the empirical 
constants employed in the ks turbel ence model- 
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In contrast lo these, not much Information Is available 

in the literature on the mathematical modelling of top injection 

CIS 1 62> 

s>^tem. However , some studies * " report advantages of bottom 

injection over top injection processes but donot address these 

C 1 

from f undamental /mat hemati cal point of view. Asai and Szekely 

were first to propose a mathematical model for the top injection 

systems. The authors formulated the flow problem by the stream 

C 1 81 

f uncti on—vor ti ci ty based method together with the 

C 1 Q1 

Pr andtl -Kol mogorov formula to calculate the turbulent 

C 171 

viscosity. Numerically predicted results gave qualitative but 

not quantitative agreement with experimental measurements reported 

earlier by Wakelin^^^^ on the CO^~water systems. Even though the 

authors reported reasonable agreement between numerical 

predictions and experimental measurements, the assumptions 

incorporated in the model were paid very little or no attention 

Cviz. , the shape of the cavity was assumed to be known and the 

vorticity considered spatially uniform, an adjustable parameter 

over the cavity surface etc.l. Furthermore, the mechanism leading 

to the transfer of momentum from the jet to the liquid bath was 

also not adequately considered and properly justified. 

C211 

Zhang et al . proposed a mathematical model for 

combined top and bottom blowing steel making process. The problem 

was formulated by statements of the turbulent transport equations 

for stream function and vorticity incorporating the k-s model. To 

describe the top injection process, the authors considered the 

momentum transfer from the impinging gas jet to the bulk liquid 

via shear stresses acting over the cavity surface CFlg. I. ID 

r 21 3 

They further related the momentum flux Ci . e. , shear stress:^ 
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across "the cavlly surface* “to the velocily gradie-nt, and the 
viscosity of the gas Csee also later D by an empirical relation : 


TCr3 


4 C Q p 

a g g ^2 

■ I -C - 3 -^ j— :> 3 


... Cl . 4:3 
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In which C , K and C are the empirical constants Cfor other 
a r 

symbols, refer to the nomenclature^ . 

The associated bottom injection phenomena was modelled by the 
quasi -single phase modelling approach. However, the plume geometry 
was not assumed to be known C rather estimated!) and gas void 
fraction accordingly was deduced from a new empirical relation: 


a Cr . z!) 


a exp 
m ^ 


2 

r 

t 3 

C H-Z+h„ 

g o 


.... c 1 . 5:3 


w3^ere a is the void fraction at the axis of symmetry and is given 
m 

as; 

a CzD * C. Q 3^^ CH-Z+h ...CI.6D 

m t g O 

Numerically predicted flow patterns were compared with 
experimental measurements and qualitative agreement between the 


two has been demonstrated. 



1.3 Scope Of The Present Work 


The brief literature review presented in section 1.2 
clearly indicates that hydrodynamic modelling of combined top and 
bottom blowing steelmaking processes has practically received 
little or no attention. Since fluid dynamics play key roles in 
determining the overall process efficiency and that adequate 
informations on flow variables etc. in high temperature 
metallurgical reactors can be usefully derived via mathematical 
model studies, it is therefore, appropriate to direct some efforts 
in this direction so that some understanding of the h>^rod>nnami cs 
of a combination blown steel making process can be developed and 
used to draw useful inferences on associated phenomena such as, 
scrap dissolution, refractory wear, slag-metal reactions etc.. 

So far, as has been pointed out already, a mathematical 

model describing the turbulent recirculatory flow in a combined 

C 21 3 

blown reactor vessel has been proposed by Zhang and coworkers 

The model predictions, however, have not been rigorously evaluated 

against equivalent experimental measurements. Experimental flow 

patterns have only been assessed against theoretical predictions 

and qualitative agreement between the two has been claimed. Since 

detailed flow measurement studies carried out on aqueous models of 

similar systems have been reported subsequently by investigators 

elsewhere, consequently the procedure adopted Cas well the 

C21 3 

assumptions incorporated in the mathematical model 3 by Zhang 

et al . needs to be re-evaluated. Furthermore, Zhang and coworkers 
have also employed a nev/ calculation procedure to model the bottom 
injection system that is much different from the procedures 
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adopted by the pr evious investigators^^ . Since numerous 
calculat-lon proc^duri&s ca.n applied t,o n^pr^serit. "the 

boll oiti/s ubmer Qe?d ga.s Inje'clion syslemi* conseque-nll y , Ib^ proc^duri^ 
adopled by Zhang el al . as well as the suggested superiorly of 
their technique above those of others require critical 
examination. This follows since no rigorous comparison between 

theory and experiments have been made to Justify the assertions 

^ C2i:> 
made 

Final ly, the investigation of Zhang and coworkers 

appears to suggest that the transfer of momentum from the 

impinging gas phase to the bulk liquid takes place only via the 

shear stress acting across the cavity surface. Consequently* the 

shearing force has been considered to be the only likely mechanism 

driving the flow field and accordingly* a hydrodynamic model based 

C 21 D 

on this concept has been proposed . However, no evidence has 
been put forward to justify such considerations. In view of the 
fact that numerous metallurgical s>'stems are driven by potential 
energy Cviz. , filling ladles, bubble stirred ladles etc.!), the 
assumption of a ** shear force driven flow " clearly warrants some 
additional considerations. It is with all these preceeding 
objectives in mind that the present research program has been 
initiated. To these ends relevant details have been summarised in 


the subsequent sections. 



Chapter II 


THEORETICAI. CONSIDERATIONS 

11.1 Assumptions In Modelling 

Combined top and bottom injection of gas into melts 
contained in industrial steel making units involve numerous processes 
that are inherently multidimensional in nature. The interactions among 
the various phases Cviz. , si ag» metal and gas2> present in the reactor 
vessel invariably lead to numerous complex phenomena, which include 
for example, waves and droplet formation at the free surface, thermal 
and particulate gradients in the melt, slag entrainment and so on. 
Consequently, to describe the gas injection induced flow in such 
multiphase and mul tidimenslonal s>^tems via an appropriate 
mathematical frame work, several idealisations have to be made and 
accordingly, numerous assumptions need to be invoked. 

In the present work, the fluid-flow model for the 
combination blown process has been derived incorporating the following 
assumpti ons , e. g. , 

Ci3 gas injection into a cylindrical vessel has been considered 
and the presence of upper slag phases in the system has been 
i gnor ed , 

Cii3 the wavy free surface at the top has been considered to be 
essentially flat, 

Ciii3 the gas is injected through a nozzle located at the base of the 
vessel as well as a lance placed above the melt. The effect of 
multiple nozzles and tuyeres have therefore, been neglected. 
Furthermore, the top as well as bottom gas Injection s>*stems within 
the cylindrical reactor have been considered to be placed along the 
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symmetry axis. Consequently, cylindrical symmetry of the system has 
been assumed and hence, variation of the flow properties along the B 
direction has been Ignored, 

Ci\0 the bubble plume formed as a result of the submerged gas 
injection has been idealised and its geometry and dimensions have been 
assumed to be known, 

CvD the shape of the cavity at the free surface resulting from the 
impinging gas jet has been assumed to be parabolic and the dimensions 
of the latter known, 

CviD the recirculating flow in the reactor has been considered to be | 

I > 

driven by the shearing forces imposed by the impinging gas jet and the 5 
buoyancy force exerted by the rising bubbles from the submerged gas ! 

injection system, ' 

i| 

CviiD the system is essentially isothermal ‘ 

and finally, 

CviiiD the Induced flow is turbulent in nature. I 

j 

A schematic of the idealised system used for the present i 

mathematical model representation is shown in Fig II. 1. 

11*2 The Gover nl ng Equat 1 ons 

Under the imposed assumption of cylindrical symmetry, the 

flow variables have been assumed to be invariant along the O 

co-ordinate and hence interms of the cylindrical polar co-ordinate 

system Cr,e and , the statement of the problem may be conveniently 

formulated via a two dimensional elliptic C©. g. »the flow is 

C 2— SD 

reel rcul at .1 ngZ) turbulent flow model Cviz. ,a set of flow equations 

in conjunction with an appropriate turbulence model D. 
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Il^E.a The flow equations 

The flow equations in essence are the turbulent 
Navi er -Stokes equations Cvii:. , an equation of continuity with two 
seperate equations for z Cthe axial D and r Cthe radial > directional 
moment um/mot i onD and are represented below : 

Equation of continuity 

du 1 dCrvD 

— ^ « Q 

r dr 


Equation of iwation in axial direction 
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Equation of motion In radial direction 
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CII . 53 


dr 


F^'s and Fg in equations CII. 33 and Cl 1. 53 are respectively the shear 
force and the buoyancy force per unit volume imposed by the top and 
bottom gas injection systems and are essentially the driving forces 
producing the flow recirculation in the combination blown reactor. 
These will be quantified Interms of appropriate expressions in the 
subsequent sectl ons . 

II. S. 2 The turbulence model 


The effective viscosity C Py3 . appearing in equations 

CII. S3 through CII. 53 follows from a turbulence model and in the present 

C 3~8D 

in\'^stigaliont the popular and the most widely used turbulence 

CISD 

model of Launder and Spalding Cviz. , the k“-£: model 3 has been 

adopted. The governing transport equations for the turbulence kinetic 
energy » k. and its dissipation rate, can be represented interms of 

the cylindrical polar coordinates as: 

Equation of turbulence kinetic energy 
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Cl 1 . 63 

in which, , the net source term can be represented as 
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Equation of dissipation rate of turbulence kinetic energy 
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CII . 103 


The auxilary relationships embodied in the above set of equations Cviz. 
equationsCII . S3 to Cl I. 93 3 are : 


^eff 


CII . 113 


p.j. ^ ^ ^ ^ 


CII. 1S3 


and finally. 


p = a + Cl-cO p^ CII. 133 

The gas volume fraction, a, appearing in equation CII. 133 is finite and 

non-zero only within the plume region Ci.e. ,0 < r < r 3. Elsewhere in 

c 


CII. 



The 


f i ve 


the dojusl n , Cr ^ , ct is zero and consequently, f>—p 

!>• 

empirical constants Cviz. ,C^ , C_, C , q- and CT". ^ embodied in 

X Cj p '^s Jc 

turbulence model were assigned the values according 

recommendations ot Launder and Spalding. These are summarised 


the k -£ 
to the 
1 n Tabl e 
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II. 3 Modelling Of The Bottom/^ubmerged gas Injection Process 

As mentioned in the proceeding sections, the bottom 

injection system in the combined blown reactor has been mathematically 

modelled via the quasi -single phase approach Cviz. , the variable 

density formulation^ assuming that the dimensions of the rising two 

phase region Cviz. .the bubble plumed is known C typically an 

empirically determined parameter 3. The buoyancy force exerted by the 

rising bubbles within the two phase region acts vertically upwards 

Cviz. .along the axial direction:) and consequently, F_, the buoyancy 

C 

force per unit vol umeC =pj^gcO has been added to the axial 

direction momentum equation Csee equation CII.33:).It is instructive to 
note here that a, gas volume fraction is finite and non-zero only 
within the two phase region. It is througli the prior presciption of 
the gas volume fraction, a, within the two phase region that the 
buoyancy force resulting from the submerged gas injection si'stem has 
been modelled in the present work. 

The gas volume fraction within the plume can be estimated 
readily applying the principle of ^lume continuity and towards this, 
numerous alternative approaches exist Csee also latter^. Typically, 
the rise velocity of the gas-liquid mixture is estimated first and 
then from the knowledge of the plume dimensions, the gas volume 
fraction, a, is conveniently deduced via either. 
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a CzD 
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. . . CII. 14D 
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if slip conditions are applied. 

In equations CII. 143 and CII. 153, r^ represents the radius of the two 
phase plume and the rise velocity of the gas -liquid mixture at any 
axial station, z. Further, U^in equation CII. 153, is the slip velocity 
and is conventionally deduced from the knowledge of average bubble 
sizes in the system. 

The density field distribution Cviz. , p in equations CII. 23, 
CII. 33 etc. 3 in the flow domain has been derived on the basis of the 
continum approach Cviz. , equation CII. 1333 and the density field 
established a priori. This is illustrated Fig II. 2. 


II. 4 Modelling Of The Top Gas Injection System 

The configuration of the impinging gas system together with 

the resultant cavity formed on the melt free surface have already been 

illustrated schematically in Fig I.l and Fig. II. 1 respectively. To 

mathematically model the flow Induced by such top gas injection 

Cl 73 j 

system, Asai and Szekely as well as Zhang and coworkers 

assumed that the geometry as well as the physical dimension of the 
cavity are known and further considered that shear stresses 



Boundary of the 
idealised two 
phase gas-liquid 


\ Gas ♦ 

^ liquid j 
^ region / 


Three types of control volumes 

Control volume A (completely inside the two phase plume) 

p - (1 ~ oC)p^ 

Control volume B (partly inside the two phase region and 
partly inside the bulk phase) 

p = v,p(i-oc) * v^Pl 

Control volume C (completely inside the bulk phase) 

P = Pl 

Fig. II. 2. The continuum approach used to define the 
density distribution in the calculation domain 




Iransmi ‘tte-d across the cavity surface essentially translates into the 
bulk phase liquid recirculation. 

To quantify the shear stresses transmitted from the 
1 mpi nql nq gas to the bulk liquid interms of an appropriate expression, 
the empirical formula suggested by Zhang and coworkers has been 

adopted in the present work and the variation of shear stress across 
the cavity surface represented according to: 


tC r 


4 C Q p 

^99 

K n 1^ d 
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3C 1 
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. . CII . 16D 


in which 1_, d, Q represent the key injection parameters, while C , 

y a 

C , K are empirical constants Cviz. , see nomenclature!). On the basis 
r 

of equation CII. 16!), the axial and the radial components of the shear 

stress at various locations on the cavity surface have been estimated 

Cby resolving rCr!) into two mutually perpendicular stress components!) 

and considering relevant area, shear layer thickness etc. , 

the correspond! I. g specific volume shear forces Cviz, , F and F D 

s f) 2, Sb , r 

estimated and added to the axial and radial direction momentumZmoti on 
equations as the representative body forces imposed by the top gas 
injection s>'^tem Cviz,, equations CII. 33 and CII. 533. 


To correlate F and F Cthe specific volume shear force 
s , r s , z 

components, see equations CII. 33 and CII. 533 with the shear stress 
TCr3, the parabolic cavity surface has been idealised interms of 
se\'eral line segments as illustrated in Fig. 11.3. Furthermore, for 
each line segment, a characteristic radial distance Cr3 has been 
chosen Cthis essenti al 1 y corresponds to the distance of the geometric 
centre of the line segment from the axis of symmetry3 and the shear 
stress estimated via equation CII. 163. Assuming the shear stress to be 



Normal face area 
to T (r), , z' 


Normol face oreo 
to T(r)r . r' 


Axial component of shear force = z TirJz 
Radial component of shear force = r^Tlrlp 


Fig. II. 3. Segmentised cavity contour and the procedure 
adopted to evaluate components of shear force / 
stress at any point on the cavity surface . 
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prcval&n't over t-he entire line segment, the two mutually perpendicular 
stress components Cviz. , along the z and the r axes respectively!! have 
been derived, multiplied by the corresponding normal area Cviz. , see 
Fig. II -33 and added directly to the axial and the radial direction 
momentum equations as the required specific volume body forces C= 
F or F in equations CII.33 etc. 3. To this end an appropriate shear 
layer thickness has been assumed Csee also later 3. 

II. 5 The Boundary Conditions 

The boundary conditions applied to the numerical solution of 
the set of partial differential equations Cviz. , equations Cl I. 13 
thrcjugh Cl I. 1033 are; 

Ci3 At the axis of the s>Tnmetry Cr=0 , O < z < L3 

v = O 

^ dk &£ 

= o ; = O and = O 

&r 9r 

Cli3 At the free surface of the liquid Cz=L, O < r < R3 

u = O 

dv dk de 

= O ; = O and = O 

3z 9z dz 

and at the side walls Cr = R, 0<z<L3 and bottom surface Cz=0, 0<r<R3 
u=0, v=:0, k=0 and £=0 

Close to the vicinity of the rigid boundaries Cviz. , the walls and the 

bottom Eurface3 >diere variations in the flow properties are steep, the 

momentum Cu and y 3 and the scalar Ck and £^ fields have consequently 

ca23 


been modelled using the routine wall functions 


Ce. g. , logarithmic 



law of wall for th«p parallel to wall velocity cojnjson^nt etc.:>. These 
have been illustrated via Fig. II. 4. 

II. 6 Numerical Solution Procedure 

Since Fg, -and ^ can be conveniently estimated from 

the procedures outlined already in sections II. 3 and 11.4, it is 
therefore readily apparent that the five principal differential 
equations Cviz., equations CII.i:), CII.E:)^ CII.4D, Cll.eZ) and CII.9D3 
together with their appropriate boundary conditions C viz. , secti on 
II. 5D can in principle be solved for the five dependent variables u, 
V, p, k and 

Thus to numerically solve the set of partial differential 

equations represented by equation CII.i:) through CII.103 together with 

the associated boundary conditions and the auxilary relations, a 

CS3D 

control volume based finite differnce procedure has been adopted 

in the present work. In this, the flow domain has been discretized 
into a number of non -over lapping control volumes, and the governing 
equations integrated around all such control volumes Cstaggered 
control volumes have been considered for the two velocity component sD 
to result into a system of algebric equations, which for any general 
variable 0, may be represented according to: 

Vp ' *e®e * Vw * 

in which, A Is the centre point coefficient Csee Fig. II. SD while A , 

P ^ 

the neighbour point C2 diniensionalD' coefficients, 

W N b 

representing the combined influences of fluid convection and 

diffusion. Further, S represents the positive component of the source 

c 



“ Free surface of liquid 

(z =L. 0 $ r $ R] 

« <5v ^ 6k ^ 6e ^ 

""'O' 6^ = ° =° 


— Axis of symmetry 
(r = 0, 2 $ L) 

^ = 0,v = 0,^ = 0 and ^ = 0 


Inside wall ““ 
[r=R, 0i$z5L) 
u =0, V = 0,k 2 0 
and € s 0 


Wall functions applied 
in the immediate 
vicinity of wall 


“Bottom wall 
(2 = 0. 0:S r^R) 
\u=0,v = 0,k = 0 and e = 0 


Fig. II. 4. Representotion of the various boundary 

conditions applied to the numerical solution 
of the governing partial differential equations 
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term. Cvlz. , of S^, etc. see equation Cl I. 3D for exampleD while 

the corresponding negative component has been embodied in the 
expression of ^p' fhe centre point coefficient. The set of algebric 
equations thus deriN-ed have been solved via the Tri -diagonal Matrix 
Algorithm, incorporating a line by line solution technique. The 
nonlinearity and coupling involved in the algebric equations C these 
follow from the governing partial differential equatlonsD have been 
taken care off by adopting to an iterative calculation procedure. 
Hence, a sequential rather than a simultaneous solution procedure has 
been embodied in the present calculation scheme. 

The pressure-veloci ty inter link age in equations Cl I. ID 

through CII.5D have been solved by an implicit finite differnce 

C 

procedure refered to as SIMPLE C Semi -Implicit Method for Pressure 

C 25D 

Linked EquationsD. For analysis of the gas-liquid region, the GALA 
technique Cthe Gas And Liquid AnalyserD was incorporated into the 
SIMPLE procedure. In this, the physical properties of the fluid 
mixture in a cell in the two phase region has been averaged on a 
volumetric basisCsee Fig. II.SD This required the conventional mass 
continuity equation to be replaced by a volume continuity equation, 
Cviz. , equation CII.IDD such that volume of fluid entering a volume 
element equals the total volume of fluid leaving the control volume. 

Since an iterative calculation procedure has been adopted, a 
convergence criterion CNormallsed ^solute Residual ^urce < O. O05Dhas 
been set on all the variables Csee also laterD , and calculations were 
carried out till WARS C Normalised Absolute Residual Sour ceD value for 
u, V, k, £ and volume continuity all fell below this stipulated value. 



II. 7 The Computer Program 


The numerical solution procedure outlined in the proceeding 

section refers to a modified version of the popular TEACH-T computer 

C 0D 

program, used extensively by Mazumdar to underscore the 

hydrodynamics of ladle gas injection operations. In the present 

C 8D 

investigation, the modified version of the TEACH-T code has been 

applied incorporating several other additional modifications 
/improvements. These for example Include, 

Ci3> adaptation to an axisymmetrlc submerged as well as to an impinging 

gas jet s'yst&m through incorporation of various force components Ce. g. , i 

F , F and F into the computational procedure:> , 

B s > z s , r 

I 

Cii3 prescription of the density field in the flow domain via the GALA ! 

1 

I 

procedure, | 

f 

! 

CiiiD incorporation of different differencing schemes Cupwind, 
exoponential etc. 2) for calculation of combined convection and 
diffusion fluxes Cviz. , the coefficients Ag, A^ etc. in equation II.17Z) 

and finally, 

C 

Civ3 a procedure Cbased on the cell porosity method 3 incorporated 
to exclude the cavity volume from the domain of calculation. 

The latter modifications will be discussed elaborately in the 
subsequent sections- The flow circuit of the computer program applied 
to the present investigation has been illustrared in Fig. II. 6. 












Chapter III 


RESULTS AND DISCUSSION 

III.l Analysis Of Sources Of Error In The Confuted Results 

To assess the validity of predicted results froir. the 

mathematical model outlined in the preceeding chapter, comparisons 
between theoretical estimates and experimental measurements are 

essential. However, prior to any comparisons it is important to 

ascertain that jxiinimuni numerical inaccuracies are involved in such 
calculated results. For a given set of equations and boundary 

conditions, the extent of assosiated error in the computed results is 
expected to depend upon several numerical parameters. These Include, 

CaD the maximum allowable residual source Ce. g. ,the convergence 
cri teri onD , 

Cb:> the distribution of grid systems in the flow domain 
and 

CcD the differencing schemes applied to estimate the convection and 
diffusion fluxes embodied in the discretization equations Ce. g. , the 
coefficients A^, Aj^, etc.>. 

It is important to recognise here that these numerical parameters 
simultaneously affect the accuracy of results as well as the 
assosiated computation times, CPU. Consequently, an optimum balance 
between accuracy and CPU time has to be considered for any effective 
error analysis. 

III. 1.1. The influence of residual source sum values on con^juted restilt 

Since an iterative calculation procedure has been adopted to 
numerically solve the governing partial differential equations, a 





convergence criterion needs to be specified. To this end, a residual 
source, defined as 

'l2 ^nb ®nb ^ ^c “ ®pl ■• ■CIII.lD 

has been estimated at each nodal point and added up over the entire 

flow domain for each Individual variable Involved Cviz. , u,v,K,£ 

etc. Obviously, when the discretization or the finite difference 

equations are completely satisfied, will be zero for all equations. 

However, since this requires an extremely large computation time to be 

achieved, a suitable convergence criterion is required to be specified 

so that the value of ^ ^ R^ ^ indicates summation over the two 

dimensional flow domain!) be less than a certain small number. However, 

smaller is this allowable residual source value, the larger would be 

the computation time. Hence, it is desirable to find out an optimum 

allowable residual source value for the present problem. 

To assess the sensitivity of computed results to the value 

of residual source, R^^, calculations were carried out for a submerged 

gas injection configuration CF_ is finite, F and F are 

a s , z s , r 

respectively zero in equations CII.3!) and CII.5DD corresponding to the 

C 273 -5 

experimental conditions of Hazumdar ei al . Ce. g. , Q=l. 667x10 

3 

m /s3 . These are shown in Table III.l for ready reference. A 

relatively large number of grid system Cviz, ,15x103 have been applied 

to encompass the flow domain CL=0. 21m and R=0. ISnD providing an 

effective grid spacing of approximately 20mm in either direction. 

C 233 

Further, the hybrid differencing scheme has been applied to 

estimate the convection and the diffusion fluxes Cviz. , the 
coefficients A , A., etc. 3. 

p N 
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Table III.l Experimental conditions of different investigators used in 
the computation of results presented in Figs. III.l through III. 5 



Experi mental 

conditions of 



Mazumdar and 

^ . ca7> 

Guthra e 

Castillejos and 

Br i macombe^ 

Lehner and co- 
Cl 

workers 

Bath height, m 

O. 21 

O. 4 

1 . 08 

Radius of the 
vessel , m. 

0. 15 

0. 25 

0. 76 

Diameter of the 
nozzl e» m 

1 . 5x1 0“^ 

6. 35x1 0"^ 


Gas flow rate» 

m /s 

Cl . 67xl0~®, 

3. 33x10“® 
and 5. Oxio"®:) 

3. 6395x1 0"'^ 

1 . 667x1 0“^ 1 

1 

1 

1 

I 

j 
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The predicted radial distribution of vertical velocity 

component at two different depths in the water model has been shown in 

Fig. III-l for three different values of ^ Cviz. ,0.001, 0.005 and 

0.025^. This appears to suggest that a reasonable value for the sum of 

the residual source is approximately 5xlo”^ and further lowering of 

this value does not lead to any Improved accuracy in the predicted 

results. However, if a sum of residual source value of approximately 

0.025 is applied in the numerical calculation procedure, significant 

accuracy in predicted results appears to be lost. Consequently, for 

all subsequent calculations reported in this thesis, an allowable sum 

of residual value equivalent to 5xlO~^ has been applied. Towards this, 

it is Important to note here that in numerous equivalent previous 

studies Csee for example reference 4> , the selection of an appropriate 

convergence criterion has neither been adequately considered nor 

explicitly specified. In addition, usage of widely different 

C 7 211) 

convergence criteria has also been noted ’ . This is particularly 

critical when comparisons between theory and experiments are to be 
made, since computed results depend considerably on such criterion. 

III. 1.2 The influence of nodal configurations on computed results 

The derivation of the discretization or the finite 
difference equations from the governing differential equations require 
numerous assumptions to be invoked Cviz. ,the mode of variation of the 
dependent variables between two nodal points etc. 3 , and consequently, 
predicted results, as one would normally anticipate, would be 
sensitive to the nodal configurations applied to the numerical solution 
scheme. However, it is well known that in the limit of a very large 
number of grid network, these assumptions and their resultant 
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0.025 
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0.001 


Dimensionless radial distance, 

Fig. III.1. Sensitivity of computations to the choice of 
residual (NARS) values illustrating the 
radial variation of the vertical component 
of flow for_ three typical NARS values 
(Q =1.67x10' V/s. R = 0.15 m, L = 0.21m, 

15 xIO grid system ). 



40 


influence on cailculaied results tend to be insignif leant and 
therefore, the numerical predictions appear to approach to the exact 
solution Cin situations where such solutions can be derivedl). 
Consequently, it is essential to establish grid independent solutions 
prior to any comparison between theory and experiments. Towards 
this, clearly an optimum nodal configuration for each individual 
problem has to be known C Indeed, a more meaningful parameter would be 
the grid spacing rather than the number of grids, since the former 
would allow one to deduce the required number of grids that would lead 
to grid independent results for any domain slzeD. 

To test these, flows in a submerged gas injection system in 

a cylindrical ladle CL=0. 21m and R=0. 16nO has been modelled 

mathematically by applying three different grid configurations viz., 

10x7, 15x10 and 17x12 grid systems respectively. ¥?hile, the hybrid 

C 233 

dif ferrencing scheme has been applied , an allowable sum of 

residual source equal to 0.005 has been considered in the calculation 
procedure. A comparison among numerically predicted radial 

distribution of \^rtical velocity component at two different depths in 
the water model is shown in Fig. III. 2. This indicates that both 15x10 
and 17x12 grid systems produce results that are practically 
equivalent. Consequently, predictions via a 15x10 grid system can be 
considered to be essentially grid independent. On the basis of this, 
grid spacing of approximately BOmm have been applied in either 
directions in all subsequent computations. It is to be noted here that 
the 10x7 grid system leads to results considerably different to those 
deduced from the two other nodal configurations. 



Vertical velocity component , (m/s) x 10 



Fig. III. 2. Sensitivity of computations to the choice of 
nodal configurations illustrating the radial 
variation of the vertical component of flow 
for three different grid systems (0=1.67x10 
mVs, R = 0.15 m, L = 0.21 m , NARS = 5x10 ) 
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III. 1.3 The influence of differencing schemes on computed results 

As mentioned in the previous chapter, the discretization 
equations embody coefficients CA^, A^^, A^ etc. D that represent the 
combined influences of fluid convection and diffusion processes. To 
estimate these coefficients Ce.g. , the convection and the diffusion 
fluxesD from the relevant velocity and diffusion coefficient values, 
several schemes can in principle be applied and these, for example, 
include : 

Cl 3 the central differencing scheme 
CBD the upwind differencing scheme 
C33 the hybrid differencing scheme 
CAy the exponential differencing scheme 
and 

C53 the power law differencing scheme 

It is to be mentioned here that except for the central differencing 

scheme, all others are known to produce physically realistic solution 

C233 

for both high and low Reynolds number flows . Consequently, the 

central differencing scheme has not been considered in the present 

analysis because of such obvious limitations. 

Since these differencing schemes incorporate various 

assumptions Cfor example, the upwind scheme ignores diffusion effects 

C 233 

for all values of Peclet number etc. 3, consequently, numerically 

predicted results are likely to be a function of the differencing 

schemes employed. Thus to test the sensitivity of computed results to 

the choice of the differencing schemes, flow calculations were 

performed via all the four schemes corresponding to the experimental 

C 273 

conditions of Mazumdar et al . Cviz. . those used in the computation 
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of Figs. III.l and III.S Csee Table and the predicted radial 

distribution of vertical velocity component at two different depths in 
the model has been compared in Fig. III. 3. This clearly indicates that 
the differencing schemes employed have practically negligible 
influence on the overall accuracy of predicted results. It is to be 
mentioned here that a 15x10 grid system together with an allowable 
residual source sum of 0.005 have been applied to derive the results 
presented in Fig. III. 3. 

Figure 111.3 further illustrates that diffusive momentum 

transfer in such systems is relatively less important* since results 
of the upwind and the rest of the differencing schemes are very 

similar. Such computed results therefore, appear to indicate that 

flows in metallurgical ladles are largely dominated by the inertial 
forces and that viscous effects Cboth molecular and turbulent!) play 
insignificant roles in affecting the overall flow pattern Induced. 
This is clearly substantiated from the results presented in Table 
III.S, which indicate that the predicted mean speed of liquid 

recirculation via all the four differencing schemes are practically 
identical C maxi mum variation ± 15^0 and eqi valent to those observed 
experimentally by Mazumdar et al . 






Fig. III.3. Sensitivity of numerically predicted vertical 

velocity component in a water model of a gas 
stirred ladle^^ illustrating the influence of diffe- 
rencing scheme employed on predicted results. 




Table III. a 
via various 
measur ements 


Theoretically predicted mean speed of liquid recirculation 

differencing schemes and their comparison with experimental 

C 273 

of Mazumdar et al . 


Differencing scheme 


Mean speed of liquid recirculationCmm/s3 


Numerical predictions Experimental values 


Upwi nd 

37. 

13 


Hybrid 

40. 

07 


Power law 

40. 

07 

34. 0 


Exponent i al 


45. 32 
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III. 2 Evaluation Of Calculation Procedures For Modelling Of The Bottom 
Gas Injection Process 

In the present study, as has been mentioned already in the 
proceeding chapters, the Eulerian quasi -single phase approach has been 
adopted to model the bottom gas injection process. Because of the 
relative ease with which such computational procedure can h>e applied, 
the quasi -single phase modelling has been widely popular and 
therefore, several versions of the procedure are available in the 
literature. In essence, the quasi -single phase modelling procedure is 
applied via a prior specification of the geometry as well as the gas 
volume fraction, a, in the plume. Since the gas volume fraction 
specified within the plume region plays key role in determining the 
overall flow field in such buoyancy driven systems, therefore, the 
various procedures available for gas volume fraction calculation are 
required to be carefully examined. To this end, any comparison among 
these calculation procedures has so far not been attempted, and 
consequently it is necessary to assess the adequacy of these by 
evaluating their relative performances, so that the most accurate 
approach can be Identified and adopted to model the submerged gas 
injection system in a combination blown reactor. 

III. 2.1 Prediction of flow parameters and average speed of bath 
recirculation via different procedures and their conqsarisons with 
experimental measurements 

To estimate the gas volume fraction, a, in the rising two 
phase plume, different group of Investigators have adopted widely 
varying concepts and based on these, four distinct types of 



computational procedures have so far emerged. These are discussed 
below In detail. 

r 3 A-y 

III.Eal.a The procedure adopted by Szekely and coworkers ' 


Szekely and 
small volume element 
differ enti al equal i on 
velocity viz. ♦ 


C3* 4!) 

coworkers through momentum balance over a 

in the plume region proposed an ordinary 
for estimation of the average plume rise 
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B Ctan XD u z 
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in which, is the average rise velocity at any axial station, z. 

The boundary conditions considered applicable to equationCIII . 1!) is 


The average rise velocity has been further related to the gas voidage 
a. according to: 


Q - rrr ^aCl -cD U 
g c s 

aCzD = 

r 

Bn f ^ U r dr 


in which, is the slip velocity and is considered to be synonymus to 

the rise velocity of a characteristic single bubble in the qui scent 
liquid^®^ and is deduced from a knowledge of the average bubble 
diameter in the system viz. , 

u s = i.oec —: > cm. 33 

Further, in the boundary condition is the superficial gas velocity 


and is deduced from : 





u 


o 


Q 


9 


n r. 


cm . 43 


Thi& void fractions, a, "thus obt^ain^d are clearly cross-sectional 
average void f'r actions. It is apparent that once the plume geometry is 
known and the void fraction distribution within the plume established 
via equations Cl II. 13 and CIII.i23» the equations of motion can be 
readily solved to yield the distribution of flow variables in the 
entire calculation domain. 


C 5 153 

The procedure adopted by Sahai and Guthrie ^ 

Sahal and Guthrie used an alternative concept to estimate 

the volume fraction distribution in the rising plume. Through detailed 

C 5 63 

hydrodynamic considerations, the authors * suggested a simple, yet 
effective formula for estimation of rise velocity in the two phase 
plume according to ; 

^1/3 ,1/4 

Q L 

U = 4.4 CIII.53 

^ 1 


and related the plume rise velocity to the gas volume fraction via 
the volume continuity principle : 


Q 


U 


av 


2 . 
nr L 
av 


ClII.a3 


in which, r refers to average radius of the plume, which is equal to 
av 

Cl/y33r^^^ Ce. g. , the cor responding radius at the mid bath depth3 . In 
equation Cl 1 1. 63 no consideration of bubble slippage has been given 
and further, the gas volume fraction has been considered to be 



spsctially uniform ov©r t-h© ©nt-lr^ volume of the two phase plume. 


I1I.2. l.c The procedure adopted by Mazumdai and Guthrie 


C28> 


C 1?0D 

Mazumdar and Guthrie modified the procedure adopted by 

C S , 6> 

Sahai and Guthrie ’ to estimate the gas volume fraction in t-he 


plume. These authors 


cze:> 


adopted the macroscopic plume model of Sahai 


and Guthrie Cvlz. , equation CIII.4I>D, and deduced the gas volume 
fraction in the plume from the following relationship : 


aCzZ) = 


Q -nr aCl-cO U 
g c s 


2 n U f "" r dr 

D 


.... cm . 5D 


Clearly, the value of is constant throughout the two phase plume 

and therefore a, the cross-sectional average void fraction is a 

function of z only. Equation CI1I.5D is different from equation 

cm. 2D since instead of U , U , a constant plume rise velocity has 

z p 

C 282) 

been used. The authors justified the use of a constant rise 

velocity rather than a variable one Cviz. , U D in equation CIII.SD, 

C 1 3D 

since experimental plume rise velocity deduced by Lehner et al . 

were found to be independent of the axial distance. The procedures 

C S B 28D 

adopted by Guthrie and coworkers * * like those of Szekely et 

also require cl prior specification of the plume geometry. 


€ 21 !) 

1II.2. l^d Zhang and coworkers formulation for estimation of gas 


void fraction in the plume 


Zhang et al . proposed a new approach to estimate the gas 

volume fraction within the plume. For the three previous approaches 
outlined earlier* the geometry of the plume has to be specified ct 
priori. In contrast, however » Zhang et al . formulation does not 



require prior specification of the two phase domain and instead, the 
gas volume fraction within the flow domain is deduced from the 
following empirical relationship: 


oCr.z:) = a exp t - — - 


J ...C1II.7D 


C CH-z+h^D 
g O 


in which, represents the void fraction at the axis of symmetry and 

is estimated from : 


2/3 , '•^L. ,1/3 -3/5 


aCz:>=C, Q^^I- 
m t 

" y 


3 


In equations Cl II. 73 and Cl I I. 83 » represents two empirical 

dimensionless constants which were assigned to values O. IS znd 9.1 
respect! vel y. 

To assess the relative accuracy of predictions via the four 
different approaches outlined abo\^ as well as their adequacy with 
reference to the mathematical modelling of submerged gas injection 
systems, flow calculations were accordingly carried out corresponding 
to the experimental conditions reported earlier by Mazumdar et al . 
Csee Table III. 13. 

Towards this, to solve equation CIII.13 and hence to deduce 

the gas volume fraction in the plume a computer program incorporating 

C293 

a fourth order Runga-Kutta method has also been developed. This is 

presented in the Appendix. The plume dimensions were estimated from the 
previously published experimental data and further a bubble diameter 
of lOmm has been considered to deduce the slip velocity. Incorporating 
the k-s turbulence model with each of the previously mentioned gas 
volume fraction calculation procedures, the distribution of the axial 
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CuD and the radial Cv:> velocity components in the flow domain were 
deduced through the solution of relevant momentum balance equa,tions 
and from these, the average or mean speed of liquid recirculation, 0, 
estimated as : 

J J* I Cu^+v^^^"^! Znr dr dz 
O O 

U = ...CIII.9D 

L R 

Bnr dr dz 



The above Integration has been carried out numerically adopting the 
Eulerian Integration procedureC viz. , assuming |Cu^+v^D^’^| to prevail 
over a control volume!) . 

Numerically predicted mean speed of liquid recirculation 

thus obtained have been directly compared with the corresponding 

C27!) 

experimental measurements of Mazumdar et al . derived for three 

different gas flow rates In a water model CL=0. 2im and R=0. ISnD of a 
gas stirred ladle system. It is to be mentioned here that the kinetic 
energy of the incoming gas has been Ignored for the flow calculation 
as a first approximation, since this is significantly lower C typically 
about 10% or less!) in comparison to the potential energy afforded by 
the gas rising through the liquid. 


TK« ccjnpariscn better: theory and experiments illustrated in 
Table III. 3. clearly indicates that the hydrod>Tiamics of the submerged 
gas injection s>'stem is adequately represented via the procedure 
suggested by Mazumdar and Guthrie . It is further apparent that 
Sahal and Guthrie’s prescription of gas voidage leads to a higher rate 
of bulk fluid recirculation. This can be anticipated since in their 
calculation procedure^®’ bubble slippage has not been considered. 

■’>1 V . •~'RARY 

1* I ■ 

'■■''r'fl 109084. 
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30. 22 
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With reference to the results shown in Table III. 2 it is also 
important to recognize here that Zhang et al ’s formulation, in 
contrast to the three other equivalent computational procedures, 
provides significantly lower estimates of mean speed of liquid 

recirculation. 

III. 2. 2 Estimates of gas volume fractions within the plume and their 

comparison with experimental measurements of Castlllejos and 

^ ^ . C30> 

Brimacombe 

Ext-ensive measurements of gas volume fractions within the 

rising two phase plume in water models of gas stirred ladles have been 

C *^02) 

recently reported by Castlllejos and Brimacombe . To assess the 
adequacy of the theoretically predicted voidage values deduced via the 
four different procedures outlined in the proceeding section, the gas 
volume fraction distribution within the rising plume were estimated 
corresponding to the experimental measurements of Castlllejos and 
Brimacombe Csee Table III. 2D. 

The experimental centreline gas voldages have been compared 

directly in Fig III. 4 against the corresponding estimates derived via 

each of the four calculation procedures assuming a bubble slip 

velocity of 0. 22m/s Cthis corresponds to a bubble diameter of lOmnO. 

There, it is clearly seen that the procedure suggested by Mazumdar and 

Guthrie*"^®^ provides the most reasonable estimates of the gas volume 

fraction distribution in the plume. Similarly, in Fig III. 5 the 

cross""S©cti ona.1 void rr&ctions de’duc^d via “Lh© f'our 

cal cula'ti on procedures have been compared wlbh equivalent, experimental 

observations of" Castillejos and Brimacombe . Once aqain adequacy or 

C 283 

the procedure suggested by Mazumdar and Guthrie is readily 



Centre line gas vodage 


• Experimental measurements of 
Castillcjos and Brimacombe^^ 

Following procedure of : 

— Mazun.dor ond Guthrie^® 

— Sohai and Suthrie^'^ 

— Szekely et al.^^^ 



Fig. Ill, 4. Theoretically estimated centre line gas voidages 
in a water model ladle and their comparison 
with equivalent experimental measurements 



Cross sectional average gas volume fraction 
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apparent. 


Referring to Figs. 111.4 and 111.5, one also finds that the 

gas volume fraction deduced via the procedure of Sahai and 

Guthrie is somewhat over estimated along almost 75% of the plume 

height. Similarly, the void fraction deduced via Szekely et al. *s 
C3,43 , 

procedure is somewhat smaller than those experimentally observed. 

In contrast to all these, Zhang et al. ’s formula^^^^ provides 

estimates of gas voldages that are significantly lower than the 

experimental values. More important, however, is their^^^^ estimates 

of centreline gas \'oldages that seem to Increase with the axial 

distance, z. which clearly contradicts the experimental findings. 

On the basis of the preceeding discussion, one would 

therefore expect that the hydrod>'namics of a submerged gas injection 

system can be simulated adequately via the calculation procedure 

C 2©D 

suggested by Mazumdar and Guthrie . While one would also expect 
that Sahai and Guthrie’s prescription of a would provide higher rates 
of bulk recirculation, equivalent procedures of Szekely and 
coworkers^^’ as well as that of Ziiang and coworkers^^^^ would 
produce smaller rates of bulk recirculation in the reactor vessel. 
Such a trend in predicted results is readily apparent in Table III. 2. 
It is instructive to note here that in all the calculation procedures 
adopted for the estimation of a, barring the procedure of Zhang and 
coworkers , radial variation of gas voidages in the plume has been 
ignored^^**^’®*®’^^ and hence the cross-sectional average void 
fraction at any axial distance can be visualised to be identical to 
the gas voidage at the axis of symmetry Cviz. ,the centre lineD. The 
results presented so far also appear to indicate that such radial 
variations can be ignored as a first approximation, so far as the 
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modelling of the bulk phase hydrodynamics is concerned. 


III. 2. 3 Prediction of rise velocity within the plume and their 
comparison with experimental measurements of Lehner et 

So far, adequacy of the calculation procedure of Mazumdar 

^ ^ , c se:> . . 

and Guthrie has been demonstrated with reference to the estimation 

of gas voidage distribution in the plume as well as the mean speed of 

bath recirculation. To assess the appropriateness of the calculation 

C 201) 

procedure further, particularly with reference to the prediction 

of rise velocity within the plume region, computations were also 

carried out corresponding to the experimental conditions of Lehner et 

C 1 3D 

al. Cthese are summarised in Table III. ID and the predicted rise 

velocity at various axial stations within the bath compared directly 

C 1 3D 

in Fig. III.B with the corresponding expserimental measurements 

Reasonable agreement between theory and experiments is once again 

readily apparent. With refernce to the results presented in Fig. 

111.6, it is to be mentioned here that a slip velocity CU^D equivalent 

to 0.4 m/s has been assumed. This corresponds to a bubble diameter of 

approximately 30mm, which is typical of such water model systems in 

the regime of the gas flow rates considered. 

On the basis of the preceedlng discussion, therefore, the 

C 20D 

calculation procedure adopted by Mazumdar and Guthrie has been 

applied to represent the bottom gas injection process in the 
combination blown system. 



Dimensionless radio! distonce, 



• Experimental measurements of 
Lehncr et al.^^ 

Present calculation applying the 

procedure of Mozumdar and Gutfirie'^ 

Fig. III.6. Theoretically predicted vertical velocity 
component within the rising plume and 
their comparison with experimental 
measurements of Lehner et al. 




Ill* 3 Modelling Of The Top Gas Injection System And Compar ison with 
Experimental Measurements 

In the preceeding sections, modelling of the bottom gas 
injection process has been considered in detail and an adequate 
mathematical representation for the gas injection induced flow 
developed. In this section, modelling approaches for the impinging gas 
Jet system currently available in the literature will be discussed and 
evaluated rigorously. Once the appropriate procedure has been 
identified, the desired mathematical model for the combination blown 
system can then be derived by considering these models for the bottom 
as well as the top gas injection s>'stems respectively. 

III. 3.1 The shearing forces acting across the cavity surface and its 
role In driving the bulk phase flow In the top gas injection system 

C 1 7 21 1> 

Earlier investigators * engaged in modelling of flows 
generated by an impinging gas Jet considered that the shearing forces 
acting across the cavity surfaces is primarily responsible for 
producing the buJ t flow reel rcul ati on. These studies, therefore, 
considered that the gas impinging over the bath surface forms a cavity 
and as the gas stream flows radially over the cavity surface, the 
outflowing gas stream drags a layer of fluid through shearing actio.is 
and this phenomenon produces the flow recirculation observed in such 
systems. In accordance with these concepts, therefore, the relevant 
shearing forces have been estimated from the gas injection parameters 
Cviz. , section II. 2. ID and embodied in the axial and the radial 
direction momentum equations as the appropriate driving body force 
term. The relationship between the gas injection parameters and the 


io 


COrr*5*5poridl ng shO’^r stra^ss Adopl^d in study Has b«M»n th^ 

same as that suggested by Zhang and coworkers . The concepts 

applied to estimate the required body forces from the empirical 

C SI D 

expression of shearing stresses ha\'e been considered briefly 

alredy in section II. 4. 

To assess the effectiveness of the above concept te. g. .that 
the flows are essentially dri'.'en via t he shear i ng forces.* as applied 
to the modelling of an axisymmetric top gas injection s>^tem, an 

experimental fluid flow study reported earlier on such s^cstems by 

C 31 D C 31 D 

Jagannathan and coworkers has been considered. These authors 

reported extensive measurements of flow parameters in a water model of 

a LD reactor via Laser Doppler Annemometer. The corresponding 

experimental conditions have been summarised in Table III. 4 for ready 

reference. The configuration of the Idealised cavity surface contour 

with reference to the finite difference grid s>'stems used for the 

present mathematical representation has been illustrated in Fig. 111. 7. 

This shows the segmentised cavity contour, < h'=- •-omp'^nerii ~ of < h*®’ 

shearing forces at any locatiorj on the cavity surface and the relevant 

normal area across which the stress components have been assumed to be 

acting. Further, the shearing forces thus estimated have been assumed 

to be prelevant over the control volumes lying immediately below the 

cavity. This implicitly implies that a shear layer thlcknisss 

equivalent to about a control volume dimension Ce. g. ,= SOuhjD has been 

considered for the present flow calculations. Incorporating adequate 

expressions for F and F in relevant equations of motion in terms 

s,z s.r 

of TCry and the appropriate grid parameters, and further considering 

Fg*0 Csee section II.2. i:>, computations were carried out corresponding 

C3i:) 

to the experimental conditions of Jagannathan et al. 



Assumed shear layer thickness 



Shearing force along r, T(r),. rAz 
Shearing force along z,r(r)2rAr 


Fig. III. 7. The shearing forces acting at a location on the 
cavity surface and the configuration of the shear 
layer thickness with reference to the finite 
difference grid system. 



Numerically predicted vertical velocity components at four 
diff€?rent depths in the model has been compared directly in Fig. 111. 8 
with the exp>er i mental observations of Jagannathaii and cov/orkers. This 
clearly indicates considerable differences between theory and 
experiments and the predicted flow is seen to be significantly smaller 
C by about a factor of 8 or so.> than the corresponding experimental 
val ues . 

This discrepancy between theory and experiments clearly 

requires some rigorous analysis and to this end the sensitivity of 

calculations to the value of 1^ C -distance of seperation between the 

lance nozzle and the bath free surface:> must first be assessed, since 

C31 D 

no value of 1^ has been specified by Jagannathan and coworkers . It 

is important to note here that calculations reported in Fig. III. 8 has 
been driven by assigning an arbitrary but likely value of 160mm to 
1^. Similarly the assumption of a shear layer thickness equivalent to 
the size of a v:ontrol vrdume may also introduce some uncertainity in 
the numer i c a 1 |“>r edi c t i ons . 

To evaluate these, sensitivity of flow predictions has been 
carried out to the choice of values of 1^ and the shear layer 
thickness. In Fig. III. 9, predicted vertical velocity component at two 
different depths in the model have been compared considering two 
arbitrary but probable values of the parameter C=150mm and 220i?:m 
respective! y3. This evidently indicates that predicted results are 
relatively insensitive to the precise value of 1^ chosen. Similarly, 
the shear layer thickness has been arbltrarly increased by a factor of 
two C shear stress assumed to be prelavent over a region that is twice 
as thick as that considered in the computation of Fig III. 83 and the 
numerical predictions repeated for the conditions of computations 




Fig. III. 8. Numerically predicted vertical velocity component at different depths and 
their comparison with experimental measurements of Jagannathan et al.^' 
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Fig. in. 9 . Sensitivity of numerical predictions to 

the values of lo illustrating the variation 
of verticol velocity corriponent at twc 
different depths in the model reactor . 




ri 


45 

corresponding "Lo FIq. Ill . 0. The result^s hsive been ill ust-ral.ed in 
Fig. III. lO where the vertical velocity components at two different 
depths have been compared for two different set of values of the shear 
layer thickness. Clearly, the assumption of a shear layer thickness in 
the range Of S0~40 mm appears to be not critical to the overall 
prediction of the flow fields. These as a consequence suggest that 
uncertain! ties assosi ated with the values of 1^ and the shear layer 
thickness are not likely to be the possible reasons for the 
discrepancy between theory and experiments illustrated in Fig. III. 8. 

111.3*2 Some additional considerations in the mathematical modelling 
of flows generated by an impinging gas jet 

The impinging gas Jet in a top blown system displaces fluid 
from the cavity region and consequently, leads to an inflntismal 
increase in the bath height. Under steady state conditions when the 
bath height assumes a constant value during the injection process, it 
is readily apparent that some energy must be continuously fed into the 
system to sustain the increased b,at h height. Towards this, a sample 
cal ciil at i on C cor r espondi ng to t he exper 1 menl al condi i i ons of 
Jagannat.han and cowor ker " :)suggests that displacement of fluid from 

the cavity region necessiates a force that is approximately equivalent 
about 1 3. BN while correspc»ndi ng cont.r ibuti on from the integrated shear 
force appears to be only 6.213x10 N. Evidently* the shearing forces 
alorie could not have produced the increased bath height and hence 
displace an eiiuivaleni amour /i of fluid frc;mi ttie* cavity regi c^ii. 
i y , it is legitimate to c-onsider 3it appiop)iaie gravity 

t or r e (.mass of the fluid disj>laced ^ acrceJ araf i on dtie to gr avi ty j in 
addi < i ori to the .shearintf fc>r*ces, as an addltiorial rhiving force 


Assumed shear layer thickness 

Below the cavity surface 

Both cbovc and below the 

cavity surface. 
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Fig. III.10. Sensitivity of numerical predictions to 
the choice of the shear layer thickness, 
illustrating the variation of vertical 
velo:ity component at two different 
depths in the model reactor. 
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producing the flow reci r culiklion In such sysieins. 

On t-htf* hasl s of thfts© concepts, at each nodal point 

immediately below the cavity surface, a representative body force has 

been evaluated C vJ , i nl er ms of mass of displacend fluid and the 

accelaration due 1.. .ji avi I y . added »o 1. he a:-d al directi on momentum 

equation i ri addi t i on to the relevant shearing force components and the 

numerical predictions carried out again corresponding to the 

C 31 y 

experimental conditions of Jagannathan and coworkers . A comparison 
between predicted and experimental vertical velocity component at four 
different depths in the model has again been illustrated in 
Fig.III.il- The i ncorp>oration of the force of impingement into the 
mathematical model is clearly seen to produce estimates of flows that 
Cln contrast to Fig. III. 83 gives excellent agreement with the reported 
experimental measurements. 

III. 4 Mathematical Modelling Of The Combination Blown Process And 
comparison Of Numerical Predictions With Experimental Measurements 

Mathematical modelling of the top as well as the bottom gas 
injection processes has been discussed extensively in the preceedlng 
chapters. Since, the appropriate mechanisms of momentum transfer from 
gas to the liquid has been identified for both top and bottom gas 
injection systems, consequently, the concepts and the procedures can 
be combined to provide an effective hydrodynamic model for the 
analysis of the combination blown steelmaking process. 

To assess the effectiveness of predictions via such an 
approach Cvlz. , coupling the models for bottom Csection III.2. 1.c3 and 
top gas Injection Csection III . 2. 23systems3 , calculations were 
performed corresponding to the experl mewntal conditions of 23iang and 
coworkers^^^^ . These exper i mental conditions are summarised in Table 



Numerically predicted vertical velocity component at different depths and 
their comparison with experimental measurements of Jagannathan et al.^' 


III. 4. In Fig. III. 12 numerically predicted flow patterns In the 
combination blown model vessel has been shown. This on comparison 

with equivalent experimental measurements of Zhang et 

C 21 3 

al . Csee Fig. III. 1 3Z) Indicate reasonable agreement between theory and 
obser vat i ons . 

In deducing the results presented in Fig. III. 12, the cavity 

volume has been artificially excluded from the domain of calculation 

by applying the *cell porosity’ concept of Moult and cowork ers^^®^. 

The principle involved has been briefly outlined in Fig. III. 14. 

The'r6*> iL is r^adiJy se*e‘n Lh^L assigning a large S value Lo ‘those 

P 

nodal points within the cavity no Flow condition can he at once 

prescribed. SindJarJy, tor the control volumes lying partly within the 

cavity arid partly within the bulk phase cell porosity factors Cf*s3 

have been estimated geometrically to deduce the effective convection 

and diffusion fluxes. Through this technique, the influence of the 

cavity surface and the cavity volume on the bulk phase hydrodynamics 

can be adequately taken into account. 

Towards these, it is instructive to note here that numerical 

C 21 3 

procedures of Zlhang and coworkers as well as that of Asai and 

C 173 

Szekely didnot consider the above and instead regarded the cavity 

as a part of the calculation domain. The cavity volume typically is 

C213 

much smaller Cfor the experimental conditions of Zhang et al . » 

cavity volume = 0.1% of the reactor volume or cavity volume = S% of 
the plume volume3 in comparison to that of the bulk and it may 
therefore appear at the first instance that inclusionT'exclusion of the 
cavity volume into/from the calculation domain may not effect the 
predicted results significantly. 








1. A completely blocked u control volume 

Sp=-10^° ond SyS 0.0 

2. A portly blocked u control volume 

— ► Aj^ = A|yj , Ag = ^ ^ *‘^2^ 

Aw'Awi ore blockage rotios 

Fig. III.U. The configuration of the cavity surface with 
reference to the finite difference grid, illust- 
rating the principles involved in the cell porosity 
method applied to exclude the cavity volume from 
the domain of calculation . 


Table III. 4 Exp>erl mental conditions of different investigators 
used in the computation of results presented in Figs. III. 7 
through III. 16. 



Experl mental 

condi ti ons of 


2Ihang and co— 

Jagannathan and co~ 


CEi:) 

workers 

. C31D 

wor k er s 

Bath hel ght , m 

0.111 

O. 324 

Radi us of the 
vessel , m 

0. 1 45 

O. 5 

Top gas flow rate. 

1 . 58x1 

8. 33x1 o“^ 

Bottom gas flow rate. 

3. 75x1 0~® 


3/s 

m 



I Finer di ameter of 

— .O 


orifice of lance 
nozzle^ m 

6x10 

2. 3x10 

Distance between 
orifice of lance 
nozzle and free 
surface* m 

0. 154 

0. 15** 

Depth of the 
bath cavity* m 

O. 0179 

O. 08 

Radius of the 
bath cavity* jn 

O. 0453 

O. 075 


Assumed 



To l^s:l lh«ps;i»» comptfl all on* w»r«i carried out again for Ihfi? 

CSl D 

experimental conditions of Zhang et al . Incorporating the cavity 

volume Into the domain of calculation. In Fig. III. 15 predicted 
ver tl cal vel oci ty component at two di f f ei eril dept hs i n the model 
vessel have been compar ed wi t h equi val ent pr edl c 1 1 ons i gnor 1 ng t he 
cavity volume. Clearly the predicted results, as one might have 
anticipated are not critical to either inclusion or exclusion of the 
cavity volume into/from the domain of calculation. It is important to 
note here that predicted vertical component of the flow deduced by 
ignoring the cavity, as seen from Fig. III. 15 is somewhat larger than 
its counterpart. This to be expected since with the exclusion of the 
cavity from the calculation scheme, the net buoyancy force acting on 
the system is accordingly also under estimated and hence this leads to 
an overall slowering of fluid motion in the system. It is also 
reasonable to anticipate here that under such conditions if the gas 
volume fraction with in the plume is increased proportionately , the 
flows deduced from either of the calculation procedure would be 
practically equivalent. 

Finally, in Fig. III. 16, the corresponding assosiated 

turbulence kinetic energy field distribution is shown. This indicates 

that higher turbulent regions C >10 ms ^ are restricted close to 

the free surface region as well as the region Immediately below the 

cavity while elsewhere, the turbulence kinetic energy is relatively 

low. Close to the bottom as well as the central portion of the reactor 

-4 2 “2 

the turbulence is seen to be extremely low C <10 m s These 

together with the predicted velocity fields may provide useful 
inferences about refractory wear, scrap melting etc. in an actual 
industrial steel making systems. 



Radial distance, mxlO 

Fig. III.15. Theoretically predicted vertical velocity 
component ot two different depths in 
the model vessel^^ illustrating the effect 
of considering /excluding the cavity 
volume on the computed results. 





Fig. III. 16. Prediction of turbulence kinetic energy (m ^ ) 
distribution in water model of o combination 
blown reactor (corresponding to the experi- 
mental conditions of Zhang et al.^’ os shown 
in Table III.4). 



Chapter IV 


CONCl.UDING REMARKS 

The following uia ior conclusions can be drawn from the 
present investigation ; 

C13 Four different procedures to model the bottom gas 
injection phenomena Cviz. ,the bulk flow, the gas volume fraction 
in the plume etc.D have been rigorously assessed and it has been 
demonstrared that the buoyancy forces Imparted by the rising 
bubbles in such systems is best represented if a constant plume 
rise velocity together with bubble slip is considered in the 
ma t-hemati cal model . 

CS:> The comparisons between the theory and the experiments 
illustrate that the shearing f or ces Imposed by the top gas 
injection are not likely to be the only mechanism producing the 
bath recircul ation in top gas injection systems and that a 
mechanism such as the force of impingement C estimated interms of 
the appropriate gravity forceD plays a crucial role in driving the 
flow field. This concept is new and is applied for the first time 
to model an impinging gas Jet induced flow system. 

C33 It has been shown that the force of impingement is many 
folds larger than the corresponding shear forces and consequently, 
an approach based on this has been adopted to model the top gas 
injection induced motion. The predicted results thus obtained are 
shown to agree reasonably well with the experimental observations 
reported in the literature. 
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C43 Numerically predicted results also indicate that 
estimates of velocity etc. in top or combination blown s>'stem are 
likeky to be relatively Insensitive to the assumed shear layer 
thickness, provided a sufficiently large number of grid points are 
applied to the numerical solui.iorj scheme. 

C53 In modelling the top gas injection induced flows, the 
cell porosity method has been applied to exclude the cavity volume 
from the calculation domain and the predicted results thus 
obtained Indicate that excluslon/lnclusion of cavity volume 
from./into the calculation domain have only marginal effect on the 
overall accuracy of predicted results. 

C63 Combining the appropriate concepts for the top and the 
bottom gas Injection procedures, a mathematical model for 
computation of flows in a combination blown steel making reactors 
has been developed and the numerical predictions evaluated against 
reported experimental measurements deduced from the water model 
studies and reasonable agreement between the two demonstrated. 



Chapter V 


RECOMMENDATIC»«S FOR FUTURE WORK 


The scope of the present work can be further expanded 
considering t.he following : 

Cl 3 The axi s>nnmetric mathematical model developed in the 
presexit study is to be applied to investigate mixing and 
mass “transfer phenomena in such systems. 

LBJ An idealised cylindrical vessel geometry has been 
considered in the development of the present model. The model must 
be made broad enough in its scope of comput.aiion so that an actual 
industria-A conf i gur Ati on C which are slightly distorted cylinder s3 
can be employed, 
and fi nally, 

C33 The axis>Tnmetric model applied in the present study can 
be modifiedy^improved to a fully three dimensional model so as to 
study theoretically the effect of multiple tuyeres or nozzle which 
are typical of industrial combination blown systems. 
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APPENDIX 



TO FIND THE SOLUTI W OF A DIFFERENTIAL EQUATION 
duCzD/dz + uCzD/z = pi wQ»*g/CaCtanClamadaD)tuCzD»z:)>«t2D BY 
RUNOA-KUITA METHOD WITH THE HELP OF GIVEN BOUI^IDARY 
CON DI TI ONS z =z Cl D , uC z =uC 0> . THE SCLUTl ON OF THE ABOVE 
DIFFERENTIAL EQUAITC^J GIVES THE VELCXIITY AS A FUNCTION 
OF z. 

CS^-ENC UNI T=23 , FI LE= ’ RUNGA. IN':) 

OPENC UNI T=22 . FI LE= ' RUNGA. OUT' 3 

IMPLICIT REAL»f8CA-H,0-Z3 

DI MENSI CM ZC 503 . UC 503 . ALFAC 503 , RC 503 

READC23.»3 C RC N3 , N=1 . 21 3 

ZC 1 3 =0. 01 842 

UC13=18. 86 

US=0. 2259 

H=0. 003335 

NMAX=21 

PI =3. 14 

Q=3. 3333E-05 

G=9. 81 

TANLAM=D. 086 
DO 20 N=1 .NMAX-1 

SI =PI *Q»*G^'C 2MC TANLAM»tUC N3 »ZC N3 3 >Hf23 -UC N3 /ZC N3 
ZCN3=ZCN3 +H/^ 

UC N3 =UC N3 +S1 *C H.."23 

S2=PI kQ»«3/C 2MC TANLAM»UC N3 WZC N3 3 »H(23 -UC N3 N3 

ZC N3 =ZC N3 +H/2 
UC N3 *UC N3 +SZitC H/23 

S3=PI 2C TANLAMitUC N3 xZC N3 3 »«t23 -UC N3 x-'ZC N3 

ZCN3=ZCN3+H 

UCN3 =UCN3 +S3WH 

S4 =PI i«Q»*G-''C 2C TANLAM»^UC N3 »ZC N3 3 3 >*)f23 -UC N3 /ZC N3 
ZCN+13=ZCN3 +H 

UC N+1 3 =UC N3 +C H/63 itC Si +BitS2 +2^S3 +S43 

CONTI NUE 

130 40 N=1 .NMAX 

A=PI »fRC N3 »«t2»«: UC N +1 3 +US3 

F=4 . OHPI »RC N3 »*t2»US)«Q 

P=2. O^PJ XRC N3 )m2>fUS 

ALFAC N3 =A.''P-C DSC5RTC A>t)t2-F3 3 .-'P 

CONTINUE 

WRI TEC 22, 4003 

FC^MATC 8X . ’ NO' , 9X , ' ZC M3 ’ , 9X . ' VELC m/s3 ’ , 8X , * RADC M3 ' 

, '6X, ' ALFACN3 '/3 
EX3 30 N=1 .NMAX 

WRITEC22,5003 N, ZC N3 , UCN3 , RC N3 , ALFACN3 
FCX?MATC5X.I4.5X.F10. 4,5X.F10. 4.5X,F10. 4,3X,F10. 6/3 
CONTINUE 
STCX= 

END 



